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o : 

i We present a theoretical analysis of inequivalent classes of interference experiments with non- 

. abelian anyons using an idealized Mach-Zender type interferometer. Because of the non-abelian 

nature of the braid group action one has to distinguish the different possibilities in which the 
experiment can be repeated, which lead to different interference patterns. We show that each setup 

■ will, after repeated measurement, lead to a situation where the two-particle (or multi-particle) state 
' gets locked into an eigenstate of some well defined operator. Also the probability to end up in such 

an eigenstate is calculated. Some representative examples are worked out in detail. 

> : 

ID . I. INTRODUCTION 

7-H ■ 

■ Physics in lower dimensions is not always simpler. A nontrivial instance in two dimensions is the possibility of 
anyons, i.e. excitations which exhibit statistics properties which transcend the conventional possibilities of bosons 
and fermions jjj. It is striking that such exotic theoretical possibilities appear to be realized in nature. With the 
experimental advances in condensed matter physics related to the fractional Quantum Hall effect, quasi-particles have 

"pH i been discovered that appear to exhibit anyonic behavior. 

There is also a more mundane interest in the two dimensional physics of anyons in the emerging field of quantum 
computation. In quantum computation, as the number of qubits grows, the problem of beating decoherence becomes 
a major threat: the larger the quantum system to be manipulated gets, the more it interacts with the (noisy) 
environment, causing the quantum state to collapse. In order to compensate the inaccuracies caused by decoherence, 
the idea of fault-tolerance and error-correction has become an important subject, and quantum codes were devised 
where single qubit-information is encoded into multiple qubits Q. If errors in the physical qubits can be corrected 
in time, the stored quantum information remains intact, or at least the decoherence time of the encoded qubit is 
increased. The progress achieved in the theory of quantum error codes however is mainly mathematical in nature. 
Current experimental realizations of this idea of qubit-encoding do not yet share the optimism of the theorists and 
are still far from being decoherence free. 

It may therefore be profitable to go back to some basic physics and search for intrinsically decoherence free quantum 
states. Promising candidates are "global quantum states" carrying quantum information that is not 'stored' locally, 
in single spins for instance, but in a robust collective property of a multi-particle system - a topological feature for 
example. Such a global characteristic of the state cannot be destroyed by local interactions, especially not by those 
with the noisy environment. As anyons exhibit such a global, topological robustness as a basic property, it seems 
anyons are especially suited for this kind of the quantum computer game, which goes by names as topological quantum 
computation || and geometrical quantum computation jj]. 

Anyons for quantum computation can be divided in two groups: the (conventional) abelian anyons which are 
manipulated by mere geometric (abelian) phases, as in ^,^], and, most promising, the so-called non-abelian anyons 
(which are strictly speaking a generalization of the abelian anyons). Non-abelian anyons may carry non-trivial internal 
degrees of freedom which can have non-abelian topological interactions which allow (in principle) for a rather clean 
manipulation of entangled states. These non-abelian anyons may serve as intrinsically fault-tolerant decoherence free 
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qubits ■ One might think of non-abelian anyons in particular media with for example broken symmetries which 

support quasi-particle excitations with anyonic properties as in model systems like the so-called discrete gauge theories 

A sometimes underestimated issue but necessary requirement for quantum computation is the ability to perform a 
measurement of the quantum state. In |^-^], the so-called fusion]] properties of the non-abelian anyons are used to 
perform a quantum measurement. But a perhaps more general applicable way to approach measurements, is to look 
at interference experiments with the quasi-particles, both for non-abelian anyons, as suggested by ||], and abelian 
anyons, |T^| . Interference experiments with non-abelian anyons are generalizations of the Aharonov-Bohm effect, and 
have been studied in the past [ fU3|]l4| ; however, these papers focussed primarily on determining the cross-section via the 
initial probability distribution, and neglected the change in the quantum state of the target particle. For topological 
quantum computation, the quantum state of the target particle, i.e., the qubit one likes to measure, obviously is 
important. 

Motivated by these considerations, we have analyzed certain types of idealized interference experiments which one 
could perform with non-abelian anyons, thereby taking into account the changes in the quantum states of the quasi- 
particles. This leads to a critical study of the non-abelian generalization of the Aharonov-Bohm effect and how it 
could be measured by a Mach-Zender type device . It turns out that there are essentially different possibilities to 
perform the repeated experiment leading to very different interference patterns. 

The paper is organized as follows, after a brief introduction to non-abelian anyons and the braid group we discuss 
in section || the general setup of experiments with the Mach-Zender interferometer. In section III we use this setup to 
analyze the different classes of anyonic experiments that could be performed and show to what kinds of interference 
patterns they may lead. We distinguish the so-called one-to-one, many-to-one and many-to-many experiments. The 
observable patterns may be very distinct, because after a sufficient number of repetitions the system will get locked 
into a particular eigenstate or subspace of some well determined operator which depends on the type of experiment. 
The proof of this locking is relegated to the appendix. The probabilities of finding certain patterns is also calculated. 
At the end of the paper we have summarized the results for the various classes in a comprehensive way in Table |. 



A. Non-abelian anyons and the braid group 



If we physically move two particles around each other and return both to their original position, this corresponds to 
a closed path in the two-particle configuration space. In three or more dimensions this closed loop will be contractible 
to a point, which means that the monodromy operator describing the net effect on the two-particle Hilbert-space 
will just be the identity operator. If we take two identical (indistinguishable) particles one may also consider the 
exchange operator defined as the square root of the monodromy operator, which consequently can only have two 
distinct eigenvalues namely plus or minus one corresponding to bosons and fermions respectively. It is well known 
that in two dimensions the situation is fundamentally different; the loop is no longer contractible so that the (unitary) 
monodromy operator can generate an arbitrary phase. Similarly the exchange operator R may yield any phase on 
a state with two identical particles which are therefore denoted as anyons. For such particles the anti-clockwise 
exchange operator 1Z and the clockwise exchange 1Z~ X have to be distinguished. 

On an algebraic level this leads to a considerable generalization because it allows also the possibility of nontrivial 
braidings between different types of particles. In short, the simple representation theory of the symmetric group 
S n (of permutations) has to be replaced by the representation theory of the braid group B n which is in principle 
an infinite group. The one dimensional unitary representations of B n are labeled by an arbitrary phase, but its 
higher dimensional unitary irreducible representations are far more complicated. In general the n-particle Hilbert 
space can be decomposed in subspaces that will carry irreducible representations of this braid group B n and if such 
representations turn out to be higher dimensional we call the particles involved, non-abelian anyons. 

The braid group B n is generated by n — 1 elements r, and their inverses. The Ti obey the Yang-Baxter equations, 
also known as the braid relations, which are: 

nn+in = n+mn+i i = l, . . . , n - 2 (1.1) 

Ti Tj = TjTi \i -j\>2 (1.2) 



1 In the discrete gauge theories |l(],|ll| the quasi-particles are allowed to fuse and their fusion-product is a gauge invariant, 
and therefore likely to be a measurable observable. However, a detailed description of the mechanism of such a measurement 
still lacks. 
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The correspondence between the Yang-Baxter equations and braids of strings of rope is best described pictorially, as 
in Fig. | 
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FIG. 1. A pictorial representation of the Yang-Baxter 
equations TiTi+iTi = t;+it;t;+i and TiTj = TjTi with > 2 

for i — 1 and j — 3; the similarity with braids of strings of 
ropes is crystal-clear. 



The way in which non-abelian anyons are usually described, is by endowing the particles with some internal degrees 
of freedom, the non trivial braid statistics of these particles can than be implemented by coupling these new internal 
degrees of freedom to a non-abelian Chern-Simons (or statistical) gauge held which mediates the appropriate type of 
topological interactions. 

As different braidings need not commute, it becomes important to be able to distinguish between them. A convenient 
way to accomplish this is to order the system in the following way: map all particles in the two-dimensional plane 
on some fixed line and number the particles (we will take the line to be horizontal and we number r from right to 
left, starting with one). If particles i and i + I pass each other on the virtual line, we have to apply the appropriate 
exchange operator: IZi if the exchange is anti-clockwise, TZ~ X if clockwise (not that after the exchange the original 
particle i now is particle i + 1 and vice versa). Every braid can be decomposed as a product of these braids of a dj acen t 
particles, IZi and 1Z~[ 1 ; also known as the braid operators, IZi and TZ~ Y clearly satisfy the braid relations (|l . f| ) , (|l.2|). 
An example is shown in Fig. 0. 




FIG. 2. An example on ordering a system with anyons and 
classifying braids. We assign a number to all six particles and 
map them on a (horizontal) line. The shown braid, which 
moves particle 3 around some of the other particles and back 
to its original position, can now be unambiguously decom- 
posed in adjacent particle exchanges; read from right to left 
this particular braid is given 



Associated with each particle is an internal space, and identical particles carry identical internal spaces. Typically, 
and this probably is the main feature of non-abelian anyons, the internal states of different particles may become 
entangled through braiding: this interaction of the non-abelian anyons is of a topological nature, independent of the 
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distance between the anyons, and it is a global property of the system, unperturbable by local interactions. The total 
system's, multi-particle, possibly entangled, internal state denoted by \ip) lies in the tensor product of the internal 
spaces of the individual particles. In such tensor products we use the same ordering of particles as introduced above, 
in that particles which are left (right) on the virtual line should appear left (right) in the tensor product. For instance, 
imagine a system with one particle of type A (with internal space V A ) and two particles of type B (with internal 
spaces V B ) which lie left of particle A. The three-particle internal state \tp) is then an element of V B ® V B (£> V A . 

The exchange caused by braid operators is made explicitly in the tensor product as well, i.e.: 1Z, TlT 1 : V B ® V A — > 
V A ®V B . We can think of 1Z and 1Z 1 as combination of some matrix and an operator a, where a only swaps the two 
vector spaces. The matrix components of 1Z and TZ~ X depend only on the types of the two particles that is operated 
upon. The square of 1Z, the monodromy operator 1Z 2 , obviously does no net swapping and can be regarded as a 
unitary matrix, of which the explicit matrix-components depend on the particles on which it operates. 

For the work we are about to present the knowledge of the braidmatrices is assumed a priori. It is clearly of 
great importance to understand the way one gets to this knowledge from an underlying algebraic structure. For 
two-dimensional physics there turns out to be an important generalization of ordinary group theory - the theory 
of quantum groups and Hopf algebra's - which provides the natural descriptions of non-abelian anyons, in that 
it generates natural representations of the braid operator R on the tensor product of its representations. For the 
Fractional Quantum Hall Effect for example, this programme has been carried out quite explicitly Jl^Jl^j . This has 
evolved into a rapidly expanding field of research on its own, we will however not make explicit use of this machinery 
in the following. 

II. THE SETUP OF EXPERIMENTS WITH THE MACH-ZENDER INTERFEROMETER 

We introduce the setup of our thought experiments using the Mach-Zender interferometer in this section. We 
will first treat an ordinary interference experiment, then move to the Aharonov-Bohm experiment, and generalize it 
further to the non-abelian anyon experiment. We will see that for non-abelian anyons, there is a nontrivial difference 
between the interference pattern and the probability distribution. 

A. Ordinary interference 

We are interested in describing interference experiments with non-abelian anyons, which basically means that we 
investigate non-abelian generalizations of the Aharonov-Bohm effect. We should realize that it is not the topological 
interactions leading to the well known (phase) factors that cause the interference, for these interactions will only 
be able to modify existing interference patterns. So, the starting point is a basic device that produces ordinary 
interference effects, as such a device, we may for example choose the Mach-Zender interferometer. We will introduce 
some terminology that we like to use in the following sections without further reference. 

The Mach-Zender interferometer, as depicted in Fig. ||, is a two dimensional device, built of two lossless beam 
splitters^ BS\ and BS2, two 100% reflecting mirrors Mi and M2, and two detectors D\ and L>2- A particle^] enters 
the apparatus at point A on the left, where it encounters a beam splitter BS\. The beam splitter will direct the 
particle along the counterclockwise path I and/or the clockwise path II. Both paths include a 100% reflecting mirror, 
M\ or M2, and both paths end up at beam splitter BS2- The particle may now emerge from the apparatus at either 
point B and be registered by detector D\ or at point C and be observed by detector D2. 

Ordinary interference is caused by the difference in relative phase of particles traveling along path I or path II, 
which is due to a difference in path length. We assume that we can adjust the apparatus in such a way that we can 
tune the difference in path length between paths I and II with high precision. We also have to include the transmission 
and reflection coefficients U^^t'^r^ of the beam splitters BSi, see Fig. f|, because these complex factors also affect 



2 We make use of the classical terminology of 'splitting a beam', but this should not distract us from the quantum nature of 
the experiment: one particle at a time may propagate through the apparatus and it is the two components of a single-particle 
wave packet that get split by the beam splitter. When probed, by a detector like Di or D2, one would find there is but one 
particle and that it followed one path, not both. 

3 We use 'particle' in a rather abstract sense; it may be any kind of particle as long as it exhibits both wave and particle 
behavior. One may, for instance, think of electrons, or photons, but in the present context also of quasi-particles with non- 
abelian braiding properties. 
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the relative phase. Since we assume the beam splitters to be perfect, its coefficients form a unitary matrix and we 
have that 

r*=ri i* = -4 \n\ + lUf = 1. (2.1) 

We will not explicitly need the lengths of paths I and II, only the relative phases that the particle will acquire when 
traversing path I or II. We will indicate the associated phase factors by e lSl and e l9u . The amplitude ^.[-Di] of the 
particle's wave packet component at detector D\ is the sum of the contributions of both paths with the appropriate 
factors: 

A[Dx] =hr' 2 e iei + nt 2 e <te ) (2.2) 

and for A[D 2 ): 

A[D 2 }=t 1 t 2 e ie * + ri r 2 e lS ". (2.3) 

The probabilities P[I?i] and P[D 2 ] that indicate the probability for the particle to be detected at Di or D 2 are the 
absolute squares of the amplitudes A[Z?i] and A[Z? 2 ]: 

P[£>i] = \hr' 2 \ 2 + |r 1 t 2 | 2 + 2Re(t 1 r^*t; e ^- e ")) , (2.4) 



P[D 2 ] = \ht' 2 \ 2 + \nr 2 \ 2 + 2Re[t 1 t' 2 r* 1 r* 2 e t( - ei - e ^j (2.5) 
= \ht' 2 \ 2 + \r ir2 \ 2 ~2Re(t l t* 2 r{r' 2 e^~ e ^) , (2.6) 



where we used ( |2.l| ) to go from ( p.q ) to (2.6). Adding eqs. ( |2.4| ) and (2.6) confirms that the total probability is of 
course preserved: 

P[D 1 ]+P[D 2 } = 1. (2.7) 

In a real experimental setup, the particles that are directed at the apparatus will not all acquire precisely the same 
phase factors e tBl or e l9u , for instance because of a slight difference in the momenta of the incoming particles. Therefore, 
we introduce a variable q and a density function p(q) to incorporate all such factors that influence the relative phase 
factors (which are due to the experimental setup). When calculating probabilities, we have to integrate over the, now 
q dependent, phase factors e* 6l ( q \ e 1011 ^: 

" e^(i)-^(i)]p(q)dq = Qe ie , (2.8) 
where Q is real and bet ween zero and one,^| and 9 is some sort of average phase difference. After integration over q, 



the probabilities of eqs. (2.4) and ( |2.6| ) then become: 

P[Di] = \hr' 2 \ 2 + \ n t 2 \ 2 + 2QRe(t 1 r' 2 r* 1 t* 2 e l9 ) , (2.9) 

P[D 2 ] = |M 2 | 2 + \ n r 2 \ 2 - 2QRe(t 1 r' 2 r* 1 t* 2 e l8 ) . (2.10) 

We will use e 1 and e 1011 in our notation for amplitudes, and Q and e 10 for probabilities, without further mentioning 
the implied g-integration. 

So far, we have only talked about probabilities for single particles to be injected and observed. Obviously, we 
want to repeat many of such single particle 'runs' and call the result the interference pattern. In other words, the 
interference pattern /[-Di] is the number of times #[D\] that a particle is injected in the apparatus and observed by 
detector D\ divided by the total number n of injected particles, and likewise for I[D 2 ]: 



4 If Q is zero no interference is observed, because the phase factors average out to zero; this is a common aspect in every 
interference experiment: the density distribution p(q) needs to be narrow enough to keep observable interference. 
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FIG. 3. The setup of the Mach-Zender interferometer for ordinary interference experiments. At 
point A a particle may be injected in the device, where it will encounter beam splitter BSi. It 
may traverse two different paths (thereby passing either mirror Mi or M2): the counterclockwise 
path I and pick up the phase factor e %Bl . or the clockwise path II and acquire e tBl1 . The particle is 
not restricted to one path exclusively: it is allowed to travel both paths at the same time. Beam 
splitter BS2 unites both paths, which now interfere with each other, and the particle emerges at 
one of two possible exits, B or C, which will be observed by a detector, D 1 or D 2 . 




FIG. 4. The coefficients U, n, t[, r[ determine the lossless beam splitter BSi (where i = 1,2); U 
indicates transmission for incoming beams from the left, r\ stands for reflection when incident from 
the left, t'i and r\ denote transmission and reflection when incident from below. 
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i[D t ] = #m, im = #m, (2 .n) 

n n 

#[D 1 ] + #[D 2 ]=n. (2.12) 

If we repeat the experiment for a large number of particles, i.e. perform many runs under the same conditions, we 
expect by definition that the observed interference patterns become equal to the probability distributions: 

lim I[Di] = P[£>i], lim I[D 2 ] = P[D 2 }. (2.13) 

n — >oo n — >oc 

On several occasions we will give examples, where for definiteness we have chosen the following values for the 
coefficients and factors that completely determine the Mach-Zender apparatus: 



T i = r 'o = 



1 , . 1 



*i = -<j=*-S> (J' = 1 > 2 ) 



y/2' J 3 V2' 



= 1, 9 = arccos(|). (2.14) 



These values were chosen such that, when substituted in eqs. (U) and fl2.10D , the probabilities P[L»i] and P[D 2 ] 
become nice, but still representative, numbers: 

P[Dx] = P[D 2 ] = ±. (2.15) 

So, in this particular example, it is by far more likely to detect the emerging particle at detector D\; only 10% of the 
incident particles is observed at detector D 2 . 

B. The Aharonov-Bohm effect: the topological phase factor 

We will now discuss the Aharonov-Bohm version of the Mach-Zender interferometer. The Aharonov-Bohm effect 
is the most logical step in going from ordinary interference experiments to interference experiments with non-abelian 
anyons, as it is both a well understood extension of ordinary interference and at the same time features a topological 
aspect. 

The Mach-Zender apparatus is extended with a physical object <E> that we put at its center, as depicted in Fig. ||, 
and we assume that this object has a topological interaction with the injected particles.^ The incoming particle now 
acquires an additional phase factor which only depends on the topology of the path (where from a topological point 
of view $ is a puncture of the two-dimensional plane) . Paths I and II arc topologically distinct and contribute the 
additional phase factors e lAl and e lXu . The amplitude Aab[-Di] of the particle's wave packet component at detector 
Di for the Aharonov-Bohm effect obviously includes these phase factors: 

Aab[A] = hr'^e^ 1 + r^2e l9lI e lAl1 , (2.16) 

and likewise for the amplitude Aab[-D2] of the component at detector D 2 . 

The probabilities Pab[Di] to observe a particle at either detector depend on the (gauge invariant) phase differ- 
ence e^ Al ~ An ) = e* A , which is the topological phase factor that a particle picks up when it circumvents $ in a 
counterclockwise way. The probabilities Pab [Di] are: 

Pab[Pi] = \hr' 2 \ 2 + \nt 2 \ 2 + 2QRe(t 1 r' 2 r* 1 t*e l9 e lX ) , (2.17) 
Pab[D 2 ] = \ht' 2 \ 2 + \ n r 2 \ 2 - 2QRe(t 1 r' 2 r* 1 t* 2 e ie e lX ) . (2.18) 



5 In case of the famous Aharonov-Bohm effect experiment where electrons are the incoming particles, $ represents a magnetic 
flux, and the electron picks up the phase factor e me * (a is some constant and e the charge of the electron) encircling the flux 
once. 
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FIG. 5. The Aharonov-Bohm effect in the setup of the 
Mach-Zender interferometer: the presence of $ contributes to 
additional topological phase factors. Following path I gives 
rise to the factor e lAl , and traversing path II yields e lA ° . The 
incoming particle still emerges at either detector D\ or D2, 
but the probabilities for either exit have been altered. The 
interference pattern now depends explicitly on the difference 
of the relative phases e lA = e i(Al ~ An) . 



As expected, the particle still emerges from the apparatus with unit probability: 

P AB [J>i] +Pab[D 2 ] = 1. 



(2.19) 



We will use the above expressions for PabA] later on, where we also want to explicitly include the dependence on 
e iA p Qr ^ii s purpose we introduce P e is\Dj\ : which is defined as: 

P e *x[A] = PA B [A]. (2.20) 

The single particle run can be repeated for many particles, which gives rise to the interference patterns Jab [Di] 
similar to those of eq. (2.11). And again, when the number n of particles grows, the interference patterns approach 



the probability distributions: 



lim J AB [A] = Pab[A] = P e "[£>i], 



(2.21) 



lim J A b[J> 2 ] = Pab[D 2 ] = PMD2]. 



(2.22) 



This is about all there is to say about the Aharonov-Bohm interference experiment apart from an exa mple that 
demo nstrat es ho w e lX changes the interference pattern. Let us substitute e lX = —1 and the values from ( [2.14 ) into 
( [2.17D and ( |2.18[ ); this yields: 



Pab[J>i] = 4j Pab[D 2 ] = ± 



(2.23) 



where the particle is now with 90% chance detected by detecto r D 2 in contrast with the values of the probabilities in 
the example of the ordinary interference experiment given by (2.15). 
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FIG. 6. The setup for interference experiments with 
non-abelian anyons using the Mach-Zender interferometer. At 
the center of the apparatus resides a non-abelian anyon, par- 
ticle A. The non-abelian anyon incident from the left is parti- 
cle B. The non-abelian anyons have a topological interaction; 
if particle B traverses the counterclockwise path this corre- 
sponds to the braid operator 1Z acting on the two-particle 
internal state, the clockwise path yields the inverse operation 
. The probabilities to emerge at detector D\ and D2 will 
explicitly depend on the expectation value of the monodromy 
operator 1Z 2 . 



C. The setup with non-abelian anyons 

The setup for - what up to now should be considered thought experiments with non-abelian anyons using the Mach- 
Zender interferometer, is very similar to the Aharonov-Bohm effect. The incoming particle will be a non-abelian anyon 
and the object $ at the center of the apparatus will be replaced by another particle which is also a non-abelian anyon. 
So, it is no longer a single particle problem (one particle scattering of an invariant external object <£•), but a two 
particle problem (one particle scattering of another one, both being non-abelian anyons). We will designate the 
incident non-abelian anyon as particle B and the non-abelian anyon fixed at the apparatus' center as particle A. The 
idea of the experiment remains the same as for the conventional interference experiments: we inject the particle B, 
in the Mach-Zender apparatus and detect it again with one of the two detectors D\ and D 2 , also see Fig. 

In calculating the amplitudes and probabilities for the non-abelian anyon experiment we have to include the two- 
particle internal state which is an element of the tensor product of the internal spaces V B and V A of particles B 
and A: 

\ip)eV B ®V A (2.24) 

. The topological phase factors in the Aharonov-Bohm effect will be replaced by braid operators acting on the state 
\ip). If particle B follows path I within the apparatus it picks up a counterclockwise braid on the two-particle internal 
state: lZ\ip); if particle B would traverse path II this would correspond to a clockwise braid: TZr 1 ^). 
Let us write down the amplitudes Ana [A] for the interference experiment with non-abelian anyons: 

AnaPi] - t ir ' 2 e^n\i>) + nhe^TZ- 1 ^), (2.25) 

Aka[D 2 ] = ht' 2 e iei n\^) + nrae^V,- 1 ^). (2.26) 

In calculating the probability distributions Pna[-Dj], we use that the braid operators are unitary and \ip) is properly 
normalized: 







Psa[Dx] = \hr' 2 \ 2 + \r 1 t 2 \ 2 + 2QRe(t 1 r' 2 r* 1 t* 2 e t9 (iP\ll 2 \^)) , 



(2.27) 



Psa[D 2 ] = \ht' 2 \ 2 + \nr 2 \ 2 - 2QRe(t 1 r' 2 r* 1 t* 2 e ie WR 2 \^)) 



(2.28) 



Compared with the probabilities Pab[Pi] for the Aharonov-Bohm effect, equations ( 2.17 ) and ( 2.18 ), the t opolo gical 
phase factor e lX has been replaced by the expectation value (ip\lZ 2 \ip) of the monodromy operator 7Z 2 in (2.27) and 

(HI). 

We can re- write the probabilities Pna[Pi] when we use that the monodromy operator is unitary and thus has 
eigenvalues of the form e tX . If we introduce the projection operators E x which project onto the eigenspace of e lA , we 
have that: 



n 2 



A E X , 



p x = (i>\E x \Tp), 



(2.29) 



irc a w>>=X> e<A > 

A 



A 



E x = 1, = !> 



(2.30) 



wh ere w e indicated the real-valued expectation value of E x by p x . With eq. (2.20) in mind, the probability distribution 
of (2.27) can now be cast in the following form: 



P NA [Pi] = \hr' 2 \ 



/ |2 



A 



pa[|*i^ 2 | 



\nt2\ 2 

I |2 , 



2QRe(^ir(,r^e 4e ^p A < 



VM 2 



2QRe(t 1 r' 2 r* 1 t* 2 e l8 e 1 ^ 



(2.31) 



This equation is nice because it gets rid of the uninteresting coefficients r^Q, etc; it states that the probability to 
observe particle B at detector D\ is a sum over probabilities of the Aharonov-Bohm version of the experiment, where 
the eigenvalues e lX of the monodromy operator 1Z 2 play the r ole of the topological phase factor. The weight p x in the 
sum is precisely the same weight as in the decomposition of ( 2.30| ) of (^}\TZ 2 \ip) into eigenvalues e lX . Of course, there 
is a similar decomposition for Pna[-D2]j so we can write for P^jDi]: 



A 



and the total probability still equals unity: 

Psa[Di] 



Pw[£ 2 ] = ]>>a^a[^] = 1. 



(2.32) 



(2.33) 



Clearly, the values of the probabilities PnaI^j] depend on the two-particle internal state \ip), via the weights p x . 
Let us gi ve an example to illustrate this. We plug the coefficients of the Mach-Zender apparat us gi ven by (2.14) into 
eq. (2.32), and assume that the monodromy operator 1Z 2 has two eigenvalues: +1,-1 (eq. (3.46) shows a specific 
example of such a matrix 7Z 2 ). In the following table, we have given Pna[-Di] for three internal states, each of which 
is characterized by the p x ; we also included the expectation value ("0|7?- 2 | tp}- 
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It is worth mentioning that one is not restricted to the Mach-Zender interferometer to make the decomposition of the 
probabilities, like in expression (2.32). This is not so strange, since any dependence on the Mach-Zender apparatus, i.e., 
the coefficients r i7 Q etc., is absent in ( 2.32| ). For instance, for double slit or plain scattering interference experiments 
with non-abelian anyons one can write down a similar decomposition. However, it is not conventional to do so; one 
usually keeps the expectation value of 7Z 2 explicit in one's expressions. 
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FIG. 7. The many-to-many experiment. To be able to mea- 
sure the initial probability distribution as the interference pat- 
terns requires that exactly the same run be repeated often. 
This means that after each observation of particle B at either 
detector both particles A and B have to be replaced by new 
ones having the same initial internal state. The old particles 
cannot be re-used because the detection of particle B at Di 
causes a change in the internal state. This way to repeat 
the runs is called the many-to-many experiment as it involves 
many particles B and many particles A. 



D. From probabilities to interference patterns 



The interference experiment with non-abelian anyons is obviously a generalization of the Aharonov-Bohm exper- 
iment, by replacing the phase factor e iA by a unitary operator 1Z 2 , which has eigenvalues of the form e lX . But the 
picture of the interference experiment with non-abelian anyons is not complete yet, and requires further analysis. We 
have to consider the question of interference patterns, which for the case of non-abelian anyons we avoided so far. 

What one wants is, to predict what one would observe in a realistic experiment. The interference pattern is 
the number of particles the detectors count if the experiment is repeated with many particles. For the ordinary 
interference and the Aharonov-Bohm experiments, is was perfectly legitimate to identify the interference patterns 
with the probability distributions, since by definition these are equal when the same experiment is repeated many 
times. Is this also true for the experiment with non-abelian anyons? The answer to this question turns out to be: 
not necessarily. The problem lies in the definition: in order to be the same one needs to repeat exactly the same 
experiment. For the non-abelian anyons experiment, this would imply that for each consecutive run, both particle A 
and B would have to be discarded, to be replaced by new particles A and B, as is also shown in Fig. [7| 

Why do we need a new particle A and the same particle B each time, why not use the old particles? The answer 
is, the two-particle internal state not only affects the probabilities to observe particle B at either detector, but is 
also altered itself by the detection, becoming a new, different state \tp'}\ It is obvious that an experiment with the 
internal state \ip) gives another result than an experiment with internal state \ip')\ so indeed, if we want to find the 
calculated probabilities as the observed interference pattern, we need a large supply of particles A and B, all with the 
same internal states Then we have that: 

lim 7£?A m [A] = lim = P NA [A], (2.34) 

n — *oo n — >oo Ti 

where 7™; m [Z)j] indicate the interference patterns for the experiment with non-abelian anyons that is repeated with 
many particles B directed to many particles A. We will refer to this type of experiment as a many-to-many (m-m) 
experiment. 
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However, we may decide to repeat the experiment in a different way, i.e., not using a new A and B each time, but 
then we will have to include the change in the internal state. The change in the two-particle internal state from 
to itself is also worth examining, and observable as well; if we want to observe the effect of this change of \ip), 
it is necessary that we repeat the experiment while keeping at least one of the two old particles A and B. We will 
discuss two such ways to repeat the experiment. First we will describe an experiment in which we keep using the 
same particle A and particle B for all consecutive runs, we call this the one-to-one experiment. The second way to 
repeat the experiment, is to keep only particle A and replace each particle B by a new one for each run; we will call 
this type of experiment many-to-one^ 

The reason that the two-particle internal state \ip) gets changed, is because the detection of particle B is a quantum 
mechanical measurement; there are two components of the two-particle wave packet, one associated with detector D\, 
t he o ther w ith D 2 , and the measurement process projects out one of these components. The amplitudes A^p^Di 



(2.25) and (2.26), describe these components: the D\ component is t\r' 2 e l lZ\i\)) + r\t 2 e "IZ" 1 ^) and the D 2 com- 
ponent is tit 2 e l<jl 1Z\'ip) + ri^e 1011 !^ 1 ^) . If particle B is detected by D\ the two-particle internal state |V) changes 
to IV'') j which is the normalized D\ component: 

|V/) = -jL (hr'^KliP) + nhe^K- 1 |V» , (2.35) 
V K 



where the real factor K is the normalization constant. There is a similar expression for the case particle B is detected 
by detector D 2 . 

A priori, we may not assume that we know the precise values of ti, Ti, r-, e 101 , e l0u in the above expression ( [2.35 ) 



it is also doubtful whether we know |V) beforehand, but even if we do, we know little of the new state |V')- The best 
we can say is that |V') will be an entangled state, i.e., i t cannot be factorized as a simple tensor product \x) 



c) £ V A , \y) £ V B ; this is because |V') is, see eq. (2.35), some linear combination of TZ and TZ 1 acting on |V), and 



for non-trivial TZ and 1Z~ 1 this action usually results in entanglement, independent of whether the initial state 
was entangled. However, as we will see in the one-to-one and the many-to-one experiment, it turns out that we can 
work with this expression for and eventually, after observing many incident particles at the detectors, we can 
even determine \ip') with precision. 



III. THE INEQUIVALENT CLASSES OF INTERFERENCE EXPERIMENTS 

In this section we want to consider the various possibilities for the interference experiments in detail. We start 
with the one-to-one experiment followed by the many-to-one experiment (we already discussed the many-to-many 
experiment in the previous section). We conclude by summarizing the different results in a Tabled . 



A. The one-to-one experiment 

The one-to-one experiment requires only one particle A and one particle B. What we will show in this subsection, 
is that after repeating the one-to-one experiment a sufficient number of times, the system will get locked into an 
eigenstate of the monodromy operator TZ 2 with a probability determined by the initial state. We inject particle B at 
the left of the Mach-Zender apparatus, with particle A at its center, and detect particle B at detector Di or D 2 , and 
we repeat this procedure using the same particles B and A. This of course requires, that we route particle B back to 
its original position - the entrance point at the left of the apparatus. We assume that we can extend the Mach-Zender 
device with some additional mirrors that take care of this, as depicted in Fig. Wf\ As we move non-abelian anyons 



One may wonder if there are more ways to repeat an interference experiment with non-abelian anyons. There are endless 
ways to do so, but these are all hybrids of the many-to-many, one-to-one and many-to-one experiments, and introduce no new 
concepts. On might be tempted to say a one-to-many experiment is missing; however, the setup is symmetrical under the 
exchange of particle A and B, and a one-to-many experiment would yield the same results as the many-to-one experiment does. 

7 These mirrors should be somewhat smart, in the sense that they may adjust their position based on whether particle B 
emerged at Di or D2, to make sure that particle B is returned at the same position independent of the detector at which it 
was observed at the previous run. 

8 Particle B might need some preparation as to become a wave packet suitable to 'split' at the beam splitter, for instance it 
may need to be accelerated. Devices, needed for such purposes, can be placed to the left of the Mach-Zender apparatus, but 
we will not bother ourselves with such devices as they do not contribute to the essence of the problem. 
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FIG. 8. The setup for the one-to-one experiment, which 
uses only one particle B and one particle A. Additional (intel- 
ligent) mirrors M have been placed around the Mach-Zender 
device to route particle B back to its original position after it 
has been detected by Di or D2 , where it will be injected again 
in the apparatus. As particle B is returned, it passes particle 
A in a counterclockwise fashion, which means the braid oper- 
ator 1Z should act on the two-particle internal state. Particle 
A remains all the time unperturbed at the center of the ap- 
paratus. With each detection of particle B, the two-particle 
internal state may change. 



around, we should be aware that this may include some action of the braid operators; as we move particle B back, 
we exchange it with particle A in a counterclockwise way, which means that an additional braid operation 1Z acts on 
the two-particle internal state 

Let us take a closer look at the one-to-one experiment, paying special attention to the internal state. We start with 
an initial two-particle internal state \ip) = IV'o); an d insert particle B in the apparatus. We t hen observe particle B at 
one of the detectors D{. The probability to emerge at detector Di is of course given by ( [2.27 ): Pna[^j] = Pna [Di] i^ D \ , 



to show explicitly the dependence on the internal state |^>o)- Since the probability for each detector will in general be 
non-zero, there are two possible outcomes for the first run: particle B is detected at D\ or at D%. If it is at detector 
D\ and we return particle B to its initial position, the new two-particle internal state \tp\) (we will use \tpj) to indicate 
the internal state after the jth run) will be: 

v K 



4= (hr'^K 2 + niae^l) | */;„>. (3.1) 
v K 



If we would find particle B at detector D2, the new state would be: 

|Vi> = -j= (ht' 2 e^Tl 2 + nr 2 e*'l) |^>, (3.2) 

where K' is a normalization constant that will in general be different from K. 

We can easily continue this for a second run. The probability to observe particle B at detector Di at the second 
run is given by PNA[-Dj]|j/>i)i where we note that \ipi), albeit implicitly, depends on the outcome of the first run. Of 
course, there are two possible outcomes for the second run, and the resulting internal state \%p2) will depend on the 
outcomes of the first and second run. We could follow this procedure for a large number of runs, say n. However, this 
means we would have to work with an internal state \ip n ) depending on the outcomes of all n runs, meaning there 
would be 2 n possible outcomes to consider! Although this would probably lead to the correct answer, the number 2™ 
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grows so incredibly fast, that this way is not advisable, and indeed we will follow another path to the answer. [] 
We are basically interested in two things after we repeated the experiment for n runs: firstly the interference 

pattern, i.e., how many times #[A] out of the total n we observed the incident particle at detector Di, and secondly, 

what we can say about the internal state \ip n ) given the interference pattern. 

The interference pattern has something to do with the probabilities Pna[A]w-}) j = 1,2, ... ,n. Since is 

different from |^>o)> which is in the general case true and independent of the outcome of the first run, we can write 

down the following inequality concerning their probability distributions: 

Pna[A]|^> T^iWA^o). (3.3) 

One can wonder if such an inequality holds also for the jth and (j + l)th run. This would mean that the probabilities 
keep fluctuating from run to run and the resulting interference pattern would become random, and thereby also 
independent on the ini tial internal state \tpo). However, it turns out that this is not so, and after a large number of 
runs the inequality of ( |3.3[ ) converges to an equality: 



Pna[A]|^„+i) = -Pna[A]|^„) for n -> oo. (3.4) 
The equality in ( |3.4| ) is automatically satisfied if \ip n +i) an d \ip n ) would differ only by a phase factor, i.e.: 

\ip n+1 ) = e ia \ip n ) for arbitrary e ia . (3.5) 



There is a class of internal states that obey eq. ( |3.5| ): the eigenstates of the monodromy operator TZ 2 . If \ip n ) would 
be an e igenstate of TZ 2 with eigenvalue e lX then, after the («.+ l)th run, the internal state |V>n+i) would be, according 
to (3.1) if the outcome of the nth run had been observation at detector D x : 

= e la \i> n ), e la = -L ( tl r' 2 e i6l e iX + r x t 2 e ie ") . (3.6) 
v K 



If at the (n + l)th run, particle B was found at detector D 2 , the new internal state would have been, using (3.2) 

1 



= e M ' |V>„), e ia ' = -j= {t x t! 2 ^e iX + r x r 2 j 6 ") . (3.7) 



The probability distributions Pna[A]|j/j„> indeed turn out to converge to a fixed value that is associated with an 
eigenvalue e %x of the monodromy operator, when the number of runs n becomes large, i.e.: 

lim Pna[A]|v-„> = \tir' 2 \ 2 + \r x t 2 \ 2 + 2QRe(t 1 r' 2 r* 1 t* 2 e te e lX ) . (3.8) 

In eq. ( |3.8| ) we still used the notation that includes the Mach-Zender apparatus' coefficients; let us write ([T^) in the 
preferred device independent notation: 

lim P NA [AW„> =P e ^[Di\. (3.9) 

n — >oo 

It is not only the probability distributions that converge, for the internal state \tp n ) also becomes fixed: as the 
eigenstate of TZ 2 that carries the eigenvalue e lA , in other words: 

lim {ip n \n 2 \ip n ) = e lX (3.10) 

n — >oo 

So, we conclude that the internal state is locked into an eigenstate of TZ 2 ; the only change subsequent runs of the 
experiment can accomplish, is an uninteresting overall phase factor e la , as follows from (3.6) and ( [3.7] ). 

But TZ 2 has in general more than one eigenvalue. What can we say about which eigenvalue e lX will come out, i.e., 
to which eigenvalue and eigenstate do P/va[A]|i/>„) and \tp n ) w hl converge? This is where the initial internal state 
= l^o) comes in; if the expectation value of TZ 2 and \tp) is given by: 



9 Actually, the first thing we did was a numerical simulation. We did not consider all 2" outcomes, but tried to simulate 
the 'real' thought experiment: guided by probabilities we followed branches in the 'tree' of all 2 n outcomes. These numerical 
simulations lead us straight to the results we are presenting here. 
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A 

then the probability that the two particle system gets locked in the eigenstate of e lX is given by p\. 

We will now consider the actual interference pattern J^[Dj] of t he o ne-to-one (o-o) experiment with non-abelian 



anyons and argue that it equals the fixed probability distribution of (3.9). Let us assume that projection, i.e., locking, 
is not achi eved until the mth run. Then the outcomes of the first m runs were governed by other probabilities than 
of e q. (|3.9| ) and the interference pattern, i.e., the collection of outcomes at the time of the mth run, will not look 
like (|3.9| ). However, after the mth run, we can repeat the experiment for an additional n runs, where n can be much 
greater than m; all these subsequent outcomes will obey ([0]) and, since n ^> m, completely dominate the interference 
pattern, rendering the first m contributions negligible: 

lim I^[Di] = P ei x[Di\. (3.12) 

n — >oo 

Not only does the monodromy operator 1Z 2 in general have multiple distinct eigenvalues, these eigenvalues can be 
degenerate as well. This degeneracy has no effect on the resulting interference pattern, which remains of the form 



of (3.12). However, the internal state does not become projected onto a particular eigenstate, but rather onto the 
eigenspace of eigenvalue e : 

lim \ip n ) = yEMt , (3.13) 

n ^°° y/(ip\E x \ip) 

where the square root in the denominator is there to normalize the state; E\ is the projection operator projecting 
onto the eigenspace of e iX . The probability p\ can also be given in terms of the projection operator E\: 

P\ = (ip\Ex\ip). (3.14) 

We have showed that eigenstates of the monodromy operator are fixed states, and said that all initial internal 
states eventually become locked. The complete (nontrivial) proof of this convergence of arbitrary initial states is 
rather lengthy, and we have relegated it to the Appendix. 

Projection is achieved after a large but finite number of runs. How do we know when we have reached such a 
projected, fixed, state? This, we determine from the interference patterns: if we can identify the interference pattern 
^na[-^»] with extreme certainty with one and only one P e ix[Dj], we know the state is projected onto the eigenspace 
of e lX . The harder it is to distinguish the P e i\[Di] for different e lX , the longer (more runs) it takes before projection 
will be realized. 

Let us conclude the one-to-one experiment with an example, with the same values ( 2.14| ) for the coefficients of the 



Mach-Zender interferometer, and a monodromy matrix 1Z 2 with eigenvalues ±1, as we used before. In the following 
table we give the interference patterns J^a [-^»] an< ^ their probability to come out of the experiment for four different 
initial internal states (which are characterized by p\ and p = -\). So, in this example, there are two possible observable 
interference patterns: either we observe nine out of ten times particle B at detector Di and one out of ten at detector 
D2 or the other way around (see the Table). Of the two situations we may end up with, we can only calculate the 
probabilities, as it should in this quantum mechanical setting. 
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FIG. 9. The many-to-one experiment: the same particle A 
remains continuously at the center of the apparatus during the 
experiment, but particle B gets replaced by a new one after 
detection. All B particles initially carry the same internal 
state. The multi-particle internal state changes with each 
observation of the current particle B at detector Di or D2 
and may become highly entangled. 



B. The many-to-one experiment: probe one particle with many others 

In the many-to-one experiment, the experiment is repeated with the same particle A but with a new particle B 
each time. This will lead to a different result, for now the system will get locked into an eigenstate of a different 
operator called U, which will be defined below. 

The setup of the Mach-Zender apparatus in this case no longer requires the extension with mirrors as needed for 
the one-to-one experiment, but it does require a large supply, i.e., a 'bag' with many identical particles B in it, on 
the left of the apparatus, see Fig. [|. We will be using the notation/convention for braids of these multiple particles 
as we introduced in section | (i.e. writing particles left in the experiment left in the tensor product, IZi interchanging 
particles i and (i + 1) anticlockwise, numbering form right to left beginning with one). The total system's initial 
internal state \tp) is an element of the tensor product of the internal spaces V B of the B particles and V A of the A 
particle: 

\ip) e V B ® ■ ■ ■ <g> V B <g> V B ® V B ® V A . (3.15) 

Once a B particle has been used, it has become uninteresting but we will not discard it completely, merely collect 
such B particles somewhere on the right of the apparatus. 

For reasons of simplicity and to give a full solution, we will subject the internal state \ip) to the following conditions: 
the initial internal state \tp) should factorize completely, i.e., the initial internal state should be unentangled, and 
furthermore, the B particles should have trivial braiding: 

ftjM=^-V)=0-<|^) = M Vz>2, (3.16) 

where Oi is the operator that only swaps particle i and (i + 1) without affecting their internal state. Applying both 
conditions to the internal state \ip), implies that l^) can be written in the following form: 

IV) = y in • ■ • y l2 y ll x a \e ln ) . . . \e i2 )\e h )\f a ), (3.17) 

where {|ej i )}is a basis for the internal space V B of the jth particle B, and {|/ Q )} a basis for V A . 

Let us start the first run of the thought experiment: we insert the first particle B in the Mach-Zender interferometer 
and o bs erve it again at either detector D\ or D 2 - The probabilities P^p\Di] for detector Di are of course given by 
( |2.31 ), (2.32), and depend on the expectation value (^ol^fl^o) °f the monodromy operator V\ which acts on particle 



A and the first particle B, with respect to the initial internal state \ipo) = \4>)- Let us assume we observe the first 
incident particle at detector D±, then the new, changed, internal state, which we indicate by l^i), is given by: 

= -4= (hr'^TZt + ntae^ftT 1 ) |^o>- (3.18) 
v K 
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Because the first particle B is now to the right of particle A, the new internal state is strictly speaking an element 
of a different vector space than that of |^o) : 



V B <g> V A - 



V 1 



(3.19) 



Furthermore, the first particle B and particle A are now entangled, so the internal state IV'i) cannot be completely 
factorized anymore: 



\^) = y*" . . .y^z m >\e t J . . .\e l2 )\f Q )\e n ), 



(3.20) 



where the components z ai1 are determined by the coefficients of ( |3.18 ) and the components x a and y % . The second 
up to the nth B particle, still have trivial braiding, which we can generalize to the jth run of the experiment by the 
statement: 



T^Mj) = IV>i)i 



i > j + 1. 



(3.21) 



So far for the first run. 

For the second run, the next particle B is injected in the Mach-Zender apparatus, and the probabilities to observe 
it at detector Di are now governed by the expectation value ( , 0ij7?||V ; l) of the current internal state and the 
monodromy operator 1Z\ that operates on particle A and the second particle B. This expectation value will in general 
be different from the expectation value of the first run: (■0i|7?.||V'i) ^ (V'ol'^ilV'o), and thus the probabilities for the 
first and second run will be different. But, just as in the case of the one-to-one experiment, the probabilities and thus 
the expectation values will converge to a final, fixed, value: 



00. 



(3.22) 



We would like to know which states obey this equality ( 3.22 ) and what the expectation value will be. These will not 
be eigenstates and eigenvalues of a m onod romy operator, as for the one-to-one experiment, but something related. 

Let us now start solving equation (3.22) if we assume t hat t he (n + l)th particle B has been observed at detector 
Di, and use that we know \4> n +i) in terms of \ip n ), as in ( |3 . 18 ) : 



(<A„ +1 |ft I 2 l+2 |V>„+i) = i( |*ir' 2 | 2 <V'»l^+i^+aftT*il^») + tir 2 r^|e i ( ei - eiI )(V'„|7e„ +1 ^ +2 ^ n+1 |V; n ) 

= (lp n \lil +1 \lp n } , 

where the normalization constant K is given by: 

K = \hr' 2 \ 2 + \nt 2 \ 2 + 2Re(t 1 r' 2 rltle l ^- e "\ij n \nl + ^ n ) 



(3.23) 



(3.24) 



These two expressions, (3.23), (3.24), may look horrid (and we can write down another pair for the case that outcome 
of the (n + l)th run is observation at detector D 2 ), but they are not that bad if we apply the braid relations to them, 
i.e., use t he fa ct that the lZ n satisfy the Yang-Baxter relations. Using (3.21) as well, two of the four expectation 
values in ( |3.23|) simplify to: 



(ip n \n n+2 n 2 n+1 Tz n l 2 \ip n ) 

K, 



(3.25) 



(Vv. \Kj r2 K 2 n+1 Ti n+ 2\i) n ) 

K, 



(3.26) 



where we introduced k to indicate the expectation value (4'n\Ti 2 l +i\'tpn} ■ Equation (3.23) will be solved if \ip n ) is such 
that: 
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(ip n \Ti n+ iTZl i+2 TZ n+1 \i>r, 



(3.27) 



mn-^nl^Ti-l^) = kk*. (3.28) 

Let us rewrite these expectation values too, by applying then Yang-Baxter equation and eq. ( |3.21 ) once more: 

{Tpn\'ll n+l 'R 2 n+2 'R n+ i\i) n ) = (-0n|7^ +1 CT„ + 2^ +1 CT„ +2 |-0„), (3.29) 
(ll>n l^n+l^n+2^n+l l^n) = (V^I^+l^n-^^+l^n-l^l^n)- (3.30) 

Although Hn+l an d <J n+2'R-ln+i <J n+2 are distinct operators, their matrix elements (i? 2 )*^ are equal: 

{^)% = {A{f a \^)\U)- (3.31) 



We will decompose the right hand side of ( 3.29 ) in a similar way in components relative to the bases {|ej)}, {|/ Q )} 
of the vector spaces V B and V A . The internal state \ip n } decomposes as: 



|Vn)=yy^ i ™- <1 |ei)|e i )|/«)|e i J...|e il ). 

Next, we let the combined operators Ti-^+i and o n +2'R- 2 l +i <J n+2 act on a basis state (ignoring the |ej 
to give: 

^ 2 l+1 a„ +2 ^ l+1 a„ +2 | ei )|e,)|/ Q ) = (i? 2 )^(i? 2 )^|e fc )|e ; )|/ 7 ). 
If we also introduce the density matrices pa and p B of particles A and _B: 

(pa)$ = z^-^zp^...^, (pb)) = y l y 3 , 
then the expectation value of ( [3.29] ) can be written as follows: 

{ip n \Tl^ +1 a n+ 2lln + xO- n+ 2\ip n ) = 

= (R 2 )iZ(R 2 )%(PB)UpB)i(pA)« 



(3.32) 
\e n ) part) 

(3.33) 
(3.34) 
(3.35) 



[(pb)UR 1 )? p ][0>b)UR 1 )%](paK 

Tr{UUp A ). 



(3.36) 



Here, in the last two steps of eq. (3.36), we introduced - as announced - an operator U that operates on V A and is 
defined as the partial trace over the B space of the monodromy operator 1Z 2 and the density matrix of particle B: 



U$ = { PB )l (Ry$, U A = Tr B (p B Tl BA ). 



(3.37) 



Let us now reconsider our original problem in terms of the U and k. We can express k as the trace over U and the 
density matrix of particle A: 



k = Tt(U Pa )- 



The condition to satisfy equality (3.23), as given by (3.27), becomes in terms of U: 

Tt(U 2 Pa ) = k, 2 . 



(3.38) 



(3.39) 



Since the trace of an operator and a density matrix can be regarded as the expectation value of that operator, we 
have that: 



(U 2 ) = {Uf 



(3.40) 



which can only be true if k is an eigenvalue of U and p A lies in the eigenspace of the eigenvalue k. 

Just as in the one-to-one experiment, the probability distributions for the incident particles converge to a fixed 
distribution associated with the eigenvalue of an operator, and the final distribution determines the interference 
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pattern. Also, the system becomes projected onto the eigenspace of the observed eigenvalue. However, the eigenvalues 
and eigenstates are no longer given by the monodromy operator 1Z 2 as in the one-to-one experiment; in the many-to- 
one (m-o) experiment we observe an eigenvalue k of the operator U and have projection onto the eigenspace of that 
eigenvalue. If we denote the probability distributions by P^°[£)j], we see that for large n they become equal to the 
interference patterns I^A^i] ■ 



lim P^piW = lim ^ A °[£i] 

n — ►oo 

\hr' 2 \ 2 + \nt 2 \ 2 +2QRe(t 1 r / 2 r* 1 t* 2 



= P*[Di], (3.41) 

where we used P re [Di] to indicate the locked probability distribution that is associated with the eigenvalue n; of 
course there are similar expressions for D 2 - As far as the projection of the internal state is concerned, there is strictly 
speaking only a projection in a subsystem of the total system's internal state, namely the subsystem describing the 
internal state of particle A, because U acts on V A only. It is the density matrix of particle A that becomes projected 
onto the associated eigenspace of k: 

E k paE k 
Tr(E K p A ) 

where E K is the projection operator of the eigenspace of k. Just as the monodromy operator TZ 2 can have multiple 
(degenerate) eigenvalues, so can U. Each eigenvalue n has a probability p K to be observed in the end; p K is determined 
by the initial density matrix p A : 

p K = Tt(E kPa ), { PA ) a p = x a x . (3.43) 

Above, we have showed that there exist states for which the probability distributions are locked, and that there 
is a chance p K to end up in each locked state. The proof of it is almost the same as that of the convergence in the 
one-to-one experiment and can be found in the Appendix. 

To solve the many-to-one experiment we introduced two conditions: the initial internal state should factorize 
completely and the B particles should have trivial braiding. We believe these conditions are too restrictive: the 
many-to-one experiment should be possible for arbitrary entangled initial states, subjected only to the condition that 
the B particles have abelian braiding: 

TZjl^k) =e iSB \^ k ) j>k, (3.44) 

where the sb in the phase factor e lSB is also known as the spin of the B particles. In this generalized many-to-one 
experiment, the interference patterns also depend explicitly on k: 

Urn I££[D i ]=P K [D i ], (3.45) 

n — >oo 

but k can no longer be regarded as the eigenvalue of the operator U ; the definition of U becomes ambiguous, because 
U is the trace of TZ 2 and the density matrix of the B particles, and in this case pb is no longer constant. We 



have not yet succe eded in proving this for the generalized many-to-one experiment, so the claim of (3.45) with the 
condition of (3.44) remains for the moment a conjecture^®\ 



Let us now give an example to illustrate the many-to-one experiment, which is borrowed from an example in [19[. 
We will take the internal space V A of the A particle to be three-dimensional and to have a basis {|1), |2), 
The internal space V B will be two-dimensional and it will have a basis |+),|— )• The monodromy operator 
1Z 2 is thus six-dimensional. Relative to the chosen bases, the monodromy operator we will use is, in a basis 
{|+)|1), |+)|2), |-)|2), |+)|3), |->|3», given by: 



10 To prove eq. ( [j.45| ) with (3.44) certainly requires some additional properties of non-abelian anyons and physical restrictions 
on the allowed internal states, especially concerning the creation and fusion of non-abelian anyons, for which the description of 
non-abelian anyons through a quantum group may be instrumental. This is beyond the scope of the present paper. 
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n 2 = 



4v3 \ 



V 



o 

1 

-1 





_ 1 
1 

4V3 



(3.46) 



/ 



This specific TZ 2 has two eigenvalues: +1, —1, which are both three-fold degenerate. Of course, it is a unitary matrix. 
We choose the incident B particles to be in the state |+), which means that their density matrix in a basis {|+), |— )} 
is: 



PB 



1 





(3.47) 



The matrix of the operator U can be simply obtained by performing the partial trace of 1Z 2 and ps, in a basis 
{|1),|2),|3)}: 



U = Ttb(TZ 2 P b) 



-5 
0" 1 
-h 



(3.48) 



So, there are two different eigenvalues k: +1 and —A. If we now plug in the values for the Mach-Zender apparatus' 
coefficients, (2.14), in the interference patterns I^]^[Di) associated with each eigenvalue, these interference patterns 
take the form given in the Table below. One outcome of this example of the many-to-one experiment is equal to the 
outcome of the ordinary interference experiment and the one-to-one experiment: nine out of ten times the incident 
particle is found at detector D\. The other possible outcome is rather different: three out of ten times is is detected 
by detector D\ and with a 70% chance the B particle is observed at detector D 2 - 





Interference patterns 


Eigenvalue k 


/n7Pi] 




+ 1 


9 
10 


1 

10 


1 


3 


7 


2 


in 


in 



The above example also demonstrates that the eigenvalues of the operator U need not have an absolute value of 
one. This would have been the case if U were a unitary operator, but U is merely a normal operator, i.e., it commutes 
with its adjoint: 

Utf = U^U, (3.49) 

which follows from the unitarity of TZ 2 and the hermicity of density matrices. Because of the definition of U in terms 
of the monodromy operator 1Z 2 , one can easily deduce that k can be written as a linear combination of eigenvalues 
e lX of 7Z 2 and that its absolute value cannot exceed one: 

K = Y J Pxe iX , 5>a = 1, => M<1. (3.50) 

A A 



C. Summary 

Let us summarize the results of the three different thought experiments we have described in this paper: the 
many-to-many, one-to-one, and many-to-one experiments. Table | gives a compact survey of the results. 

The topological interactions of non-abelian anyons are generally agreed upon to be observable in interference 
experiments. We have demonstrated that such interference experiments are generalizations of interference experiments 
of the Aharonov-Bohm type. However, whereas in the Aharonov-Bohm versions, the probability distributions are 
identified with the interference patterns, this cannot be readily carried over to experiments involving non-abelian 
anyons. There it is crucial to explicitly describe the way in which each run of an experiment is repeated. 
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TABLE I. Summary of the results for the three different thought interference experiments with non-abelian anyons we 
discussed. Notice that the results are independent of the exact setup of the Mach-Zender interferometer and may be applied 
to other setups as well. The one-to-one and many-to-one experiment lock the system into an eigenpattern or eigenvalue 
associated with an operator (either the monodromy operator TZ 2 , or U: the partial trace over TZ 2 and the density matrix of 
particle B) acting on the internal state, and project the internal state onto the eigenspace. The many-to-many experiment 
merely determines the expectation value of the monodromy operator 1Z 2 , and ignores any change in the internal state due 
to the measurements. Although interference patterns cannot be identified with probability distributions, there is still some 
conservation of probability: the average of the outcomes of all three different schemes is the same. For ordinary interference 
experiments or Aharonov-Bohm experiments there is no need to distinguish between the many-to-many, the one-to-one and the 
many-to-one schemes, as all schemes will yield the same results in these cases, the difference becomes important for non-abelian 
anyons only. 





many-to-many 


type of interference experiment 
one-to-one 


many-to-one a 


The measured expectation/eigen value, 

{lp n \Tln+l \tpn) for 71 ^ OO 


im 2 \4>) 




K 


The diagonalised operator 


not relevant 


TZ 2 


U = Ttb{1Z 2 P b) 


The observed interference pattern, 

Aa[A] = lmin^oo fW [A] 175.2^1^) 




A- [A] 


A A] 


Change of internal state 


not relevant 




. , . ApaA 




pA Tta(Apa) 


Probability of eigenvalue/pattern 


1 


Px = (V|A|V> 


Pk = TT A (E K p A ) 


Average of all outcomes 




E A PAe lA = (m 2 \i>) 


E K p K ft = {ip\lZ 2 \ip) 


Aharonov-Bohm pattern Ab A*] 
with ft 2 = e iA l 


A* [A] 


A-a [A] 


A- A] 


The plain interference pattern I[Di] 
with TZ 2 = 1 


A A] 


A [A] 


A [A] 



a Recall that to solve the many-to-one experiment we posed some restrictions on the internal state of the B particles; the 
operator U can only be denned if the initial internal state is completely unentangled. 



In the many-to-many experiment, each new run is conducted with fresh particles A and B carrying the initial 
internal state \ip). In this case, by definition, the interference pattern reproduces the probability distribution, and is 
characterized by the expectation value (tp\7Z 2 \^) of the monodromy operator TZ 2 . 

The scheme of the one-to-one experiment uses the same two particles A and B over and over again. The resulting 
interference pattern now does not depend on (ip\TZ 2 \ip) , but rather on an eigenvalue e lX of the unitary monodromy 
operator TZ 2 . Because there is a clear analogy between the eigenvalues e lX and the associated interference patterns 
P e i\ [Di] we could use the term eigenpattern to indicate the eigenvalue-pattern P e i\ [A]- So, in the one-to-one exper- 
iment we observe an eigenpattern, or in other words we measure an eigenvalue. The measurement of the eigenvalue 
projects the two-particle internal state onto the eigenspace of e . The probability to observe a specific eigenpattern 
is determined by py. the fraction of the internal state \ip) that carries eigenvalue eA 

The many-to-one experiment describes the way one intuitively (and conventionally) envisages the repetition of an 
interference experiment with non-abelian anyons: fix particle A inside the interferometer and direct many identically 
prepared particles B at it. However, the many-to-onc turns out mathematically to be the hardest of all three schemes 
to analyze, as we need to deal with the total system's internal state each run. If we restrict ourselves to initially 
unentangled internal states, we can identify the resulting eigenpattern A [A] with the eigenvalue n of an operator U 
that acts on the internal space of particle A] the operator U explicitly depends on the internal state of particle B, as 
U is the partial trace over TZ 2 and ps- The possible eigenvalues n of U thus also depend on pb- The projection of the 
system onto the eigenspace of k effectively only affects the A particle: it is its density matrix pa that gets projected. 
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Interference experiments, including thought experiments, obviously need a device that can create interference; the 
topological interaction of the non-abelian anyons does not cause interference, it only alters existing interference. We 
chose to use the Mach-Zender interferometer and explicitly used its determining coefficients ti, r' 2 , Q etc., in many 
equations and our examples. In the end however, we switched to a notation that is independent of the exact setup, and 
only depends on eigenvalues and eigenpatterns. The results in Table | are thus not restricted to experiments with the 
Mach-Zender interferometer, but apply to interference experiments in general, for instance double slit experiments. 

Although the three schemes yield different results, the probability distribution for the first incident particle is the 
same for each scheme, and we expect to see some conservation of this initial probability. This is achieved by taking 
the average of all outcomes. When averaged, the many-to-many, one-to-one, and many-to-one experiments all return 
the same result which reflects the equal probability for the first run of either experiment. 

Being a generalization, the non-abelian anyon experiments should of course cover the Aharonov-Bohm effect and 
the ordinary interference experiments as well. If we let the monodromy operator act on a trivial internal space, 
i.e., 1Z 2 — e lA l or just TZ 2 = e lA , we retrieve the Aharonov-Bohm effect: the interference pattern for all three 
schemes will be the same P e tx[Di], and the question of projection becomes irrelevant on the trivial internal space. 
If we let 1Z 2 act as unity, we recover the ordinary interference experiment, and the three different ways to repeat 
the experiment obviously yield the same interference patterns. The difference between interference patterns and 
probability distributions becomes non-trivial only in the case of non-abelian anyons. 



APPENDIX: PROOF OF LOCKING 



Here in the Appendix, we will give the proof that in the one-to-one and many-to-one experiments, all initial internal 
states converge to a locked state. We will begin from the point of view of the one-to-one experiment, and we will 
see later on that the proof for the many-to-one experiment is essentially the same. We will repeat some of the more 
important equations from sections |lj] and [II, and we will state what needs to be proven from the physical point 
of view. From there on, the proof is purely mathematical, and we will change notation as to clearly discriminate 
physics from mathematics. The pivot of proof is to construct a new probability distribution function, which has a 
clear behavior as a function on n, and shows that \ip n ) becomes projected for large n with the correct probability. 

So, focusing on the one-to-one experiment now, we will first determine an expression for (ipnlE^lipn) as a function 
of the outcomes D^. (j S {1, . 
state after the rth r un i s 
becomes, eqs. (3.1),(|3^): 



, n}, ij £ {1, 2}) of the performed runs of the experiment. If the two-particle internal 
then based on the outcome of the (j + l)th run, D, j+1 , the new internal state IV'j+i) 



I'/'., 



+1/ 



(Al) 



when the outcome is observation at detector D\ ; if the incident particle emerges at detector D 2 , the new state is given 
by: 



1 



3 + 1 1 



(ht'^TZ 2 + r 1 r 2 e Wu t)\i> J ). 



(A2) 



The dependence on the outcome Di j+1 is not explicit yet in the above expressions. However, we can rewrite the parts 
(. . . TZ 2 + . . . 1) in terms of E\ and e lA , if we use that (eqs. fl2.29D,Q2.30|)): 



for then: 



TZ 2 =J2^E X , 1 = ^£a, 

A A 

hr'^lZ 2 + r x t 2 e iBll t = ^{hr'^ e lX + ri t 2 e i6n )E x 

A 



Here, we introduced /a [A] for the sake of abbreviation, but also the dependence on both the eigenvalues e lX and 
the outcomes Di becomes more clear. If we take the absolute square of f\[Di] (and implicitly integrate over q) we 
recognize the Aharonov-Bohm probability distribution P e i\ [Di]: 
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\t 1 r' 2 e ie *e a +r 1 t 2 e ie *y = |/ A [A]| = P e «*[I>i]. 

If we write |^> n ) in terms of \ipo), according to eqs. ( |Al[ ). ( |A2| ) and by substituting f\[Di], we obtain a product of i 
similar factors: 



Most of the projectors E\ j cancel each other, if we apply to them the projector algebra: 

E\E^ = Sx^Ex- 

Now, we are ready to express the expectation value (V , n|A*|V'rc.) with explicit dependence on the outcomes {A,} 
of the n runs, in the following way: 

W>n|A#„) = -ip^lAxl-.-Pei^Aj^ol^l^o)- (A3) 
A 

The normalization constant K is simply given by: 

K = J2 P^ [AJ ■ • • [D in ] (MExHo) (A4) 

A 

= p n (D h ,D i2 ,...,D in ), (A5) 

but this K turns out to have a physical meaning too: it is equal to the probability p n ( Ai , A 2 , . . . , A n ) to observe 
the particular sequence of n outcomes at the detectors Ai , A 2 , . . . , A„ . With this probability in mind, proving 
that the state locks means we need to show that for large values of n the following is true for a certain eigenvalue e IM : 

(4>n\E^n) -» 1, 

W£aW-»0, VA^/x, (A6) 

or equivalently: 

(</>JAl<M » W>n|Wn), VA ^ jU. (A7) 

The probability^ for this particular eigenvalue e v to come out is determined by the initial state \tpo): 

p^ = (MEM- (A8) 

We will now assume that there are only two distinct eigenvalues e lAl and e tX2 ; later on, we will give a generalization 
for cases with more than two eigenvalues. As the following part of the proof is purely mathematical, it is convenient 
to use an abstract notation which does not directly relate to the thought-experimental setup of the interference 
experiments. Let us introduce the functions A(i), B(i), p(ii, ■ ■ . , i n ) and the constants a 2 , ft 2 , by making the following 
identifications: 

A(i) = P e i\ 1 [Di] (A9) 
B(i) = P e » 2 [Di] (A10) 
p(ii, ...,i n )= p n (Ai,- • A„) (AH) 
a 2 = {^\ExM (A12) 
f3 2 = (ip \Ex 2 \ipo). (A13) 

All are non-negative, and can be easily shown to obey the following relations: 

5>« = 1, = 1, a 2 +ft 2 = l, 

i i 

p(»i, . . . , in) = a 2 A{i{) . . . A{i n ) + (3 2 B(h) . . . B(i n ), 
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^2 Pih, •■■,««) = 1. 

Next, and this is the pivot of the proof, we define a new probability distribution function P n (z) of a real- valued 
variable z, which in turn is defined through the relation: 

{i>n\ExMn) ^ ' 

i.e., the probability that after the nth run the expectation value of E\ 1 is equal to that of E\ 2 times the exponent of 
z. The formal definition of P n (z) involves a Dirac delta-function: 

Pn (z)= £ P (i 1 ,...,i w )g(*-in ( Ais ) 

As a probability distribution P n {z) should be normalized, as it is: 

P n (z)dz = 1. (A16) 



Convergence, i.e., locking, is achieved when either (ip n \E Xl \ip n ) > {ip n \E\ 2 \ip n ) or (^ n |.E\ 2 |i/> n ) > {ip n \E\ t \tp n ), which 
means that P n (z) needs to become zero near z = 0, because the area under P n {z) around z = stands for the 
probability that (ip n \E Xl \ip n ) is approximately equal to (ip n \E\ 2 \ip n ) . 
If we use that: 

(V-nl-EAilVO , a 2 A{h) . . . A{i n ) 
m -; ; — r = m — - + In — ; — - ; — - , 

(ip n \Ex 2 \ip n ) (3 2 B(h)...B(i n y 

then P n {z) can be conveniently written as the sum of two functions that are scaled by a 2 and (3 2 : 

P n (z) = a 2 P^(z-ln^) + /? 2 P,f(z-ln^). (A17) 
These two functions Pn( z ) an d Pn( z ) are independent of a 2 and /3 2 : 

P A (z) = J2 A(h)...A(i n )5(z- In (A18) 



P n s (z) = £ B(h)...B(i n )6(z-ln ^Hfej )' (A19) 



and can both be regarded as probability distribution functions, since: 



OO POO 

■A/„\J. _ 1 / nB 



P n A (z)cb = l, / i^(«)dz = l. (A20) 

Next, we will calculate the expectation values of z and z 1 with respect to P^{ z ) an d P^{z). 
The expectation value (z)^ of z for P±{z) is: 



/OO 
-OO 



B(i) 

= m A , (A21) 

where ttia > because ^ B(i) and:[^] 



A(i) unequal to B(i) means P e i\ 1 [Di] 7^ P e t\ 2 [Di], which is usually true unless Q is zero, in which case there is no observable 
interference, and thus the whole experiment would be uninteresting. 
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a In — + b In — = (a — b) (In a — In 6) > 0, a ^= b. 
b a 

The expectation value of z 2 is given by: 

(*)f = j%^)z 2 dz = £A(z)ln 2 ||j 

= 4 + m^, (A22) 
which also defines s\. We can define m B and s 2 ^ in a similar way: 

J-oa j #W 

= -m B < 0, (A23) 



<z 2 



J-oo i a 



'(<) 

= s 2 B + m B . (A24) 



If we now calculate the expectation values (z)^, (z 2 )^, etc., for arbitrary n by plugging in eqs. ( |A18D and ( |A1S| ) 
we find these depend on n in the following way: 

/oo 
P 7 f(z)zdz = n{z)f = nm A - (A25) 
-OO 

(z)% = nm B , (z 2 )f = ns 2 B + n 2 m 2 B (A27) 

The expectation values of z and z 2 tell a great deal over the general form of Pn(z) an d P^(z), by considering the 
mean \i = (z) and the variance a 2 = (z 2 ) — (z) 2 of P„ (z) and P 7 f(z): 

fJ-An = (z)n = nm A , (A28) 
<T 2 An = (z 2 )*-((z)A) 2 =ns 2 A , (A29) 

(j, Bn = -nm B , a Bn = ns 2 B . (A30) 

Both the mean fi and the variance a 2 of Pn(z) an d Pn(z) scale linearly with n. The standard-with cr, however, is 
given by the square root of the variance a 2 , and thus scales with s/n. This means that if we compare P^(z) and 
P^ Q (z), the mean of Pf^ (z) is a hundred times greater than that of P^(z), while the width has only increased by a 
factor ten. The total area under both Pf(z) and P^ Q (z) remains equal to one. 

So, P~-(z) and P r f(z) are both probability distribution functions, of which the peak 'runs away' from the origin 
when n increases, and these peaks run harder than they broaden. This behavior is in fact sufficient as proof. For let 



us look at P n (z), eq. ( A17 ), which is after all the probability distribution function we are physically interested in, for 
increasing n. The area under P n (z) for large values of z indicates the probability that after n runs (ipnlEx^ipn) = 1; 
the area under P n (z) for large negative z is the probability that (ifj n \E\ 2 \ip n ) = 1. By now, we know that P„(z) is 
completely dominated by the P£ (z)-part for large z and by the P^ (z)-part for large negative z, because of the way 
the means and widths scale with n. 

Nevertheless, let us specify "large z" by choosing a fixed z' > OQand calculate some z'-dependent areas under 
P„(z) when n goes to infinity: 



2 To be precise, if z > z' , this means that {ip n \E\ 1 \il} n } > — ■ So, if z = 25, then e z » 10 10 , and thus {il) n \E\ x \ip n ) > 
- 10 -10 : the state is pretty much projected, and we can regard z = 25 as large. 
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lim / P n (z)dz = 0, (A31) 

/>OC 

lim / PJz) dz = a 2 , (A32) 

,wo ° J z > 

lim / P n {z) dz = [3 2 . (A33) 

Clearly, the results of these integrals do not depend on our arbitrary choice of z' . And so, the projection is proved, 
i.e., at the end of the experiment the system is locked with complete certainty, including the correct probabilities for 
each possible outcome, since a 2 = p\ ± and /3 2 = p\ 2 . 

The proof for more than two eigenvalues is a simple generalization. For the case of three eigenvalues, P n depends 
on two variables z and w, where P n [z, w) then stands for the probability that at the same time both: 



(V>n|-E A3 |Vn)' (^n\E\ 3 \^ n ) 



(A34) 



P n (z, u>) can now be written as a sum of three other scaled probability distribution functions P^(z, tv), P„(z, w), and 
(because of the third eigenvalue) P^(z,w), which all 'run away' from the origin with increasing n. This procedure 
can be easily generalized for an arbitrary, but finite, number of eigenvalues. 

For the many-to-one experiment, the recipe of the proof is simply: replace the eigenvalues and projector operators 
by their appropriate counterparts, i.e., replace each e lX by k, each P ei \[ Di i by P K [Di], and each E\ by E K . Although 
\ip n ) in the many-to-one experiment is very different from \tp n ) in the one-to-one experiment, i.e., an element of a very 
different vectorspace, and the same applies to E\ and E K , we can nevertheless in the many-to-one experiment come 
up with an expression similar to ( |A3| ): 

(il>n\EM = ^P K [D il ]...P K [D in ](i> \E K \il) ). (A35) 

A 

We arrive at ( |A35| ) by determining (ipj+i\E K \ipj+i) in terms of (ipj\E K \ipj) , by doing some rewriting of 

(i> j+1 \K 2 j+2 \iJ j+1 ) = « (^i+il^l^+i) (A36) 



From eq. ( |A3q ) forward, the proof of the one-to-one and many-to-one experiments is exactly the same. 

The functions Pn( z ) an d P 7 f(z) have other interesting properties as well, for instance P„ (z) = e z P^(z), but also 
p 2 = p \ ° p \ wh ile P 2 + P\ o Pi. Further more the sums over i, where i could be 1 or 2 corresponding to detector 
D\ or D2, can easily be generalized to more detectors, i.e., sums where i can take on more values, and can also be 
replaced by integrals for setups with a 'continuous spectrum of detectors'. These aspects we will not discuss here, as 
they lie outside the current scope; see (ljj for some more details. 
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